< 
OO 
[tel 


CN 1127844 


(19) 中 华人 民 共 和 国 国 家 知识 产权 局 


ap (12) 发 明 专 利 申请 


(21) 申请 号 202110090133.2 
(22) 申请 日 2021.01.22 


(71) 申请 人 南京 理工 大 学 
地 址 210094 江苏 省 南京 市 孝 陵 卫 200 号 
(72) 发 明 人 KE AR TEE 4AN 
MEE ” 侯 晓 明 
(74) 专利 代理 机 构 南京 理工 大 学 专利 中 心 
32203 


代理 人 汪清 


(51) Int.Cl. 
GO6F 30/23 (2020.01) 
GO6T 17/00 (2006.01) 


(54) 发 明 名 称 
一 种 弹 炮 耦合 有 限 元 参数 化 建 模 方法 
(57) 摘要 


本 发 明 提 出 了 一 种 弹 炮 耦合 有 限 元 参数 化 
建 模 方法 ,涉及 机 械 结构 参数 化 建 模 技术 领 ] 
将 线 膛 身 管 弹 炮 耦合 有 限 元 建 模 过 程 通过 
Python 语言 进行 描述 ,使 用 脚本 建立 弹 炮 耦合 有 
限 元 模型 ,使 得 建立 弹 炮 耦合 有 限 元 模型 更 加 方 
便 、 快 捷 ; 主要 包括 :首先 建立 身 管 及 弹丸 二 维 草 
图 ,将 二 维 草 图 实现 全 约束 并 定义 参数 化 变量 ， 
根据 二 维 草图 自动 生成 三 维 模型 ,并 根据 缠 角 变 
化 规律 采用 Python 语 言 对 身 管 节点 进行 偏 移 , 采 
用 Python 描述 建立 弹丸 与 身 管 之 间 的 耦合 关系 。 
本 发 明 采 用 编程 技术 实现 了 弹 炮 耦合 有 限 元 模 
型 的 自动 更 改 与 数值 模型 的 自动 更 新 ,从 而 提高 
了 弹 炮 耦合 有 限 元 模型 的 建 模 速度 ,缩短 了 设计 
周期 , 极 大 提高 了 工作 效率 。 
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1 .一 种 弹 炮 耦合 有 限 元 参数 化 建 模 方法 ,包括 以 下 步骤 ; 
步骤 1 ,根据 身 管 实际 结构 尺寸 ,确定 结构 参数 变量 ,利用 Python 语言 定义 结构 参数 化 
变量 ， pre edic 
又 2 , Ti ds ORAL EET SE os Z8 T4] , TECH TR zo CIE Hj E vr ELE 2 Ep E HOA, — HERE, 
me 三 维 实体 模型 ,从 建 模 过 程 中 软件 后 台 目 动 生成 的 rpy 文 件 中 提取 出 相应 程序 段 ， 
得 到 里 管 及 弹丸 三 维 实体 建 模 程 序 ; 
h 使 用 Python 语言 控制 直 膛线 身 管 及 弹丸 的 有 限 元 单元 网 格 的 自动 化 划分 ; 
步骤 4, 导出 直 膛 线 身 管 有 限 元 inp 文 件 , 利 用 Python 语言 读 取 身 管 单 元 节点 并 进行 单 
元 节点 偏 和 ,生成 具有 螺旋 膛线 的 身 管 有 限 元 模型 inp 文 件 ; 
步骤 5 ,利用 Python 语言 恋 取 步 骤 4 生 成 的 inp 文 件 并 将 模型 导入 到 有 限 元 软件 ; 
通过 在 有 限 元 软件 ABAQUS 中 对 身 管 及 弹丸 模型 进行 材料 属性 定义 、 装 配 定义 、 
分 析 步 定义 、 接 触 约 束 定 义 、 载 和 荷 及 边界 条 件 加 载 操作 ,将 操作 过 程 中 软件 后 台 自 动 生成 的 
rpy 文 件 中 的 程序 语句 提取 出 来 ,得 到 身 管 和 弹丸 耦合 关系 施加 的 Python 程序 语句 ; 
步骤 7 ,在 第 6 步 装配 体 模型 的 基础 上 ,编写 静态 与 动态 分 析 作 业 自 动 创建 程序 语句 ; 
步骤 8 ,将 建 模 过 程 程序 进行 集成 整合 ,建立 弹 炮 耦合 有 限 元 模型 脚本 ,通过 脚本 运行 
HA ,建立 弹 炮 耦合 有 限 元 模型 。 
2. 根 据 权 利 要求 1 所 述 的 弹 炮 耦合 有 限 元 参数 化 建 模 方 法 ,其 特征 在 于 ,所 述 步骤 2 中 
的 直 膛 线 身 管 炮 口 端面 膛线 草图 采用 全 约束 , 且 膛 线 边 界线 之 间 通 过 角度 进行 约束 。 
3. 根 据 权 利 要 求 1 所 述 的 弹 炮 耦 合 有 限 元 参数 化 建 模 方法 ,其 特征 在 于 ,步骤 4 进行 单 
元 节点 偏 移 ,具体 包括 以 下 步 又: 
4.1) 定义 身 管 装配 体 ,在 装配 体 模块 定义 身 管 装配 体 模 型 ; 
4.2) 导出 直 膛 线 身 管 inp 文 件 , 将 身 管 装配 体 模 型 的 Jop-SG.inp 文 件 导出 ; 
4.3) 计算 身 管 膛线 缠 角 沿 炮 膛 轴线 的 变化 规律 , 渐 速 段 缠 角 正切 值 变 化 规律 为 y a), 
等 齐 段 缠 角 正切 值 变 化 规律 为 y', (x) : 


tan a, 一 tan a l tan a; 
y (x)= ! 2| x+ ? (0<x<1) 
l tan 2, — tan ar, 


y; (x) 2 tana, (I « x X L) 


其 中 ,x 为 膛线 上 菜 一 点 与 阴线 起 始点 之 间 的 距离 ,膛线 初始 缠 角 为 a@, 炮 口 缠 角 为 a ， 
膛线 总 长 为 L, 渐 速 段 长 度 为 ]; 

4.4) 按照 缠 角 变化 规律 进行 身 管 节点 偏 移 。 

4. 根 据 权 利 要 求 3 所 述 的 弹 炮 耦合 有 限 元 参数 化 建 模 方 法 ,其 特征 在 于 , 线 膛 身 管 有 限 
元 单元 节点 偏 移 过 程 为 : 

Python 程序 首先 从 炮 尾 中 心 点 开始 读 取 吴 管 单元 层 坐 标 和 节点 坐标 值 , 当 坐标 值 xv 小 
于 阴线 起 点 与 药 室 底 端面 之 间 的 距离 h 时 ,单元 节点 坐标 不 变 , 当 单元 层 x, 坐 标 值 大 于 h 且 
小 于 ht+1l 时 ,根据 缠 角 变化 规律 公式 得 到 身 管 渐 速 段 日 元 节点 偏 移 后 的 坐标 (x0,y z^) ; 当 
单元 层 x 坐标 值 大 于 h+1l 时 ,进行 等 齐 段 身 管 节 点 的 偏 移 , 得 到 等 齐 段 身 管 单元 节 点 偏 移 后 
的 坐标 xo. y^ z^) , 随 着 单元 识别 层 编号 i 的 增加 ,单元 节点 按照 缠 角 变化 规律 沿 着 身 管 轴 
线 方向 依次 进行 偏 移 , 直 到 单元 层 x, 坐 标 值 大 于 身 管 长 度 L. 时 ,号 管 单元 节点 偏 移 结束 ,得 
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到 节点 偶 移 后 的 混合 膛线 身 管 有 限 元 模型 Job-SG_NEW.inp 文 件 。 

5. 根 据 权 利 要 求 1 所 述 的 弹 炮 耦合 有 限 元 参数 化 建 模 方 法 ,其 特征 在 于 ,步骤 8 建立 弹 
炮 耦 合 有 限 元 模型 ,具体 过 程 为 :首先 设置 好 变量 的 初始 值 ,然后 打开 ABAQUS 初 始 化 界面 ， 
点 开 RunScript 脚 本 运行 指令 ,选中 弹 炮 耦 合 参数 化 脚本 程序 py 并 点 击 运行 ,ABAQUS 软 件 将 
按照 程序 代码 ,依次 完成 建立 实体 模型 、 划 分 网 格 、 定 义 并 赋予 材料 属性 、 设 置 接触 与 约束 
关系 、 施 加 载荷 与 边界 条 件 。 

6. 根 据 权 利 要求 1 所 述 的 弹 炮 耦合 有 限 元 参数 化 建 模 方法 ,其 特征 在 于 ,所 述 装配 后 的 
弹丸 与 映 管 之 间 设 置 接触 属性 ,接触 行为 采用 罚 函 数 法 进行 描述 。 

7. 根 据 权 利 要求 1 所 述 的 弹 炮 耦合 有 限 元 参数 化 建 模 方 法 ,其 特征 在 于 ,所 述 身 管 有 限 
元 单元 节 扣 根据 膛线 缠 角 变化 规律 进行 偏 移 。 

8 .根据 权利 要 求 1 所 述 的 弹 炮 耦合 有 限 元 参数 化 建 模 方法 ,其 特征 在 于 ,所 述 静 态 与 动 
态 分 析 作 业 之 间 通 过 预定 义 场 进行 联系 与 信息 导入 。 
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一 种 弹 炮 耦合 有 限 元 参数 化 建 模 方法 
技术 领域 


[0001] “本 发 明 属 于 结构 参数 化 建 模 技 术 领 域 , 特 别 是 一 种 弹 炮 耦 合 有 限 元 参数 化 建 模 技 
术 领 域 。 


背景 技术 

[0002] ”弹丸 发 射 是 一 个 极其 复杂 的 动力 学 过 程 。 弹 丸 在 膀 内 运动 期 间 与 炮 膛 之 间 相 互 作 
用 , 弹 带 材料 在 高 温 高 压 下 发 生 高 度 非 线 性 的 力学 行为 ,为 了 研究 弹丸 膛 内 时 期 的 动力 学 
规律 ,需要 建立 弹丸 与 身 管 相 互 作 用 的 耦合 模型 。 采 用 有 限 元 法 建立 弹 炮 耦合 模型 对 火炮 
进行 发 射 动 力学 仿真 分 析 是 研究 火炮 系统 射击 过 程 内 在 规律 的 一 种 重要 手段 。 传 统 建 模 方 
法 非常 耗 时 , 且 建 模 精 度 较 低 ,为 了 能 够 深入 研究 弹 炮 相互 作用 的 力学 规律 ,往往 需要 建立 
大 量 的 仿真 模型 进行 动力 学 仿真 分 析 。 目 前 的 建 模 方法 难以 保证 快速 建立 出 较 高 精度 的 弹 
炮 耦 合 有 限 元 模型 ,研究 建立 具有 较 高 精度 的 参数 化 有 限 元 模型 对 火炮 发 射 动 力学 仿真 具 
[0003] ”由 于 线 膛 身 管 内 膛 参 数 众多 ,弹丸 与 身 管 耦合 作用 极其 复杂 ,加 之 弹丸 与 身 管 之 
间 复 杂 的 力学 行为 容易 导致 模型 计算 不 收敛 ,这 使 得 建立 出 具有 较 高 精度 的 弹 炮 耦合 有 限 
元 模型 变 得 十 分 困难 。 现 有 技术 解决 弹 炮 耦合 有 限 元 建 模 的 手段 大 多 都 是 先 在 三 维 建 模 软 
件 中 分 别 建立 身 管 及 弹丸 三 维 实体 模 型 ,然后 将 模型 导入 网 格 前 处 理 软件 中 进行 网 格 的 划 
分 ,最 后 将 划分 好 的 吴 管 网 格 模型 导入 有 限 元 软件 中 ,最 后 进行 弹丸 与 身 管 的 耦合 关系 的 
施加 ,从 而 建立 出 弹 炮 耦合 有 限 元 模型 。 采 用 以 上 建 模 方法 存在 的 缺点 在 于 :1) 由 于 软件 兼 
容 性 问题 ,三维 实体 模型 导入 网 格 划分 软件 后 会 经 常 发 生 精度 问题 ,导致 有 限 元 网 格 划分 
失败 。2) 由 于 弹丸 与 届 管 模型 厢 合 关系 需要 进行 手动 逐步 操作 ,导致 建 模 过 程 操 作 步 又 繁 
3t FERT «3) 对 于 不 同 结构 参数 的 弹 炮 耦 合 有 限 元 模型 ,每 当 结 构 参 数 改变 时 ,都 需 重新 进 
行 所 有 建 模 步 又 的 操作 ,导致 建 模 过 程 自 动 化 程度 极 低 。 

[0004] “建立 线 膛 身 管 的 精确 有 限 元 模型 是 研究 弹 炮 耦 合 问题 的 关键 。 现 有 技术 中 ,能 够 
较 好 地 解决 复杂 空间 曲线 特征 有 限 元 建 模 的 方法 为 节点 偏 移 法 ,如 以 公布 的 发 明 专 利 
CN107679341A 一 种 身 管 结构 有 限 元 参数 化 建 模 方法 即 采 用 节点 偏 移 法 对 线 膛 身 管 进行 参 
数 化 建 模 , 但 该 方法 只 能 建立 映 管 的 有 限 元 参数 化 模型 ,未 能 解决 弹丸 与 身 管 耦合 的 参数 
化 建 模 问 题 。 


发 明 内 容 

[0005] ”本 发 明 的 目的 在 于 提供 一 种 弹 炮 耦 合 有 限 元 参数 化 建 模 方法 ,以 解决 传统 弹 炮 耦 
合 有 限 元 建 模 效 率 低 的 问题 .通过 有 限 元 软件 二 次 开发 可 以 实现 弹 炮 耦合 有 限 元 复杂 模型 
的 自动 ,快速 高效 建立 ,无需 借助 第 三 方 软件 ,只 需 在 ABAQUS 软 件 中 运行 脚本 程序 , 即 可 实 
现 弹 炮 耦 合 有 限 元 模型 的 自动 化 建立 。 

[0006] ”实现 本 发 明 目 的 的 技术 解决 方案 为 : 

[0007] ”步骤 1 ,根据 身 管 实际 结构 尺寸 ,确定 结构 参数 变量 ,利用 Python 语 言 定义 结构 参 
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数 化 变量 ,编写 有 限 元 软件 建 模 初 始 化 程序 ; 

[0008] ”步骤 2 ,根据 弹丸 及 号 管 实际 结构 ,在 有 限 元 软件 中 建立 直 腾 线 身 管 及 弹丸 二 维 草 
图 ,进而 建立 三 维 实体 模型 ,从 建 模 过 程 中 软件 后 台 自 动 生 成 的 rpy 文 件 中 提取 出 相应 程序 
段 ,得 到 身 管 及 弹丸 三 维 实体 建 模 程序 ; 

[0009] ”步骤 3 ,使 用 Python 语言 控制 直 膛 线 吴 管 及 弹丸 的 有 限 元 单元 网 格 的 自动 化 划分 ; 
[0010] ”步骤 4, 导出 直 膛线 吴 管 有 限 元 inp 文 件 , 利 用 Python 语言 读 取 身 管 单元 节点 并 进 
行 单元 节点 偶 移 ,生成 具有 螺旋 膛线 的 身 管 有 限 元 模型 inp 文 件 ， 

[0011] 步骤 5 ,利用 Python 语言 读 取 步 骤 4 生 成 的 inp 文 件 并 将 模型 导入 到 有 限 元 软件 ; 
[0012] — 2 96 ,通过 在 有 限 元 软件 ABAQUS 中 对 吴 管 及 弹丸 模型 进行 材料 属性 定义 、 装 配 定 
义 、 分 析 步 定义 、 接 触 约 束 定义 、 载 荷 及 边界 条 件 加 载 操作 ,将 操作 过 程 中 软件 后 台 自 动 生 
成 的 rpy 文 件 中 的 程序 语句 提取 出 来 ,得 到 身 管 和 弹丸 耦合 关系 施加 的 Python 程序 语句 ; 
[0013] “步骤 7 ,在 第 6 步 装配 体 模型 的 基础 上 ,编写 静态 与 动态 分 析 作 业 自 动 创 建 程序 语 
句 ; 

[0014] “步骤 8 ,将 建 模 过 程 程序 进行 集成 整合 ,建立 弹 炮 耦 合 有 限 元 模型 脚本 ,通过 脚本 
运行 指令 ,建立 弹 炮 耦合 有 限 元 模型 。 

[0015] “本 发 明 与 现 有 技术 相 比 ,其 显著 优点 是 : 

[0016] (本 发 明 极 大 提高 了 弹 炮 耦合 有 限 元 的 建 模 的 效率 ,将 弹丸 与 身 管 有 限 元 耦合 
建 模 过 程 采 用 程序 语言 编写 出 来 ,相关 人 员 在 建立 弹 炮 耦合 有 限 元 模型 时 ,只 需 将 程序 化 
脚本 在 有 限 元 软件 ABAQUS 中 运行 后 , 即 可 快速 建立 出 有 限 元 模型 , 极 大 简化 了 建 模 过 程 , 使 
得 模型 生成 速度 非常 快速 。 

[0017] ”本 发 明 实 现 了 身 管内 腾 结 构 参数 的 参数 化 设置 ,对 于 不 同 结构 参数 的 建 模 ,只 需 
更 改 参 数 化 变量 值 , 即 可 快速 建立 出 有 限 元 模型 ,同时 与 手动 逐步 操作 建 模 过 程 相 比较 ,由 
于 建 模 过 程 使 用 程序 化 操作 流程 ,从 而 降低 了 有 限 元 模型 出 错 的 概率 。 

[0018] — (2) 本 发 明 提 高 了 弹 炮 耦合 有 限 元 模型 的 精度 , 建 模 过 程 由 于 避免 了 不 同 软件 之 
间 的 交互 转换 ,从 而 使 得 有 限 元 模型 的 精度 大 大 提高 。 
[0019] — (3) 本 发 明 建 立 模型 的 种 类 广泛 ,该 方法 不 仅 适用 于 混合 膛线 里 管 弹 炮 厢 合 建 模 ， 
对 于 渐 速 膛线 和 等 齐 膛线 身 管 弹 炮 契合 建 模 同样 适用 。 


— 


附 图 说 明 

[0020] ”图 1 是 本 发 明 的 参数 化 建 模 方法 具体 实施 操作 流程 。 
[0021] ”图 2 为 号 管内 膛 结 构 示 意图 (单位 :mm 。 
[0022] ”图 3 为 号 管 外 形 结 构 示 意图 (单位 :mm 。 
[0023] ”图 4 为 膛线 横 剖 面 结构 示意 图 (单位 :mm) 。 
[0024] ”图 5 为 膛线 展开 曲线 示意 图 。 

[0025] ”图 6 为 膛线 截面 二 维 草图 (单位 :mm) 。 
[0026] ”图 7 为 直 简 圆柱 身 管 有 限 元 模型 示意 图 。 
[0027] ”图 8 为 弹丸 结构 示意 图 。 

[0028] ”图 9 是 弹丸 有 限 元 网 格 模型 示意 图 。 
[0029] ”图 10 为 喘 管 节点 偏 移 程序 逻辑 图 。 
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[0030] ”图 11 是 节点 偏 移 后 的 线 膛 身 管 有 限 元 模型 示意 图 。 

[0031] ”图 12 为 膛线 部 分 有 限 元 网 格 模型 示意 图 。 

[0032] ”图 13 为 有 限 元 软件 ABAQUS 初 始 化 界面 。 

[0033] ”图 14 为 有 限 元 软件 ABAQUS 脚 本 运行 界面 。 

[0034] ”图 15 为 弹丸 与 身 管 耦合 的 有 限 元 模型 示意 图 。 

具体 实施 方式 

[0035] ”下面 结 合 说 明 书 附 图 对 本 发 明 做 进一步 的 介绍 。 

[0036] ”结合 图 1 至 图 15 对 本 发 明 内 容 进 行 具体 介绍 ,运用 该 参数 化 建 模 方法 建立 弹 炮 夺 
合 有 限 元 模型 ,具体 操作 步骤 如 下 : 

[0037] ”步骤 1,Python 语 言 定 义 参数 化 变量 

[0038] ”首先 对 建 模 初 始 化 程序 进行 定义 ,语句 如 下 : 

[0039] from abaqus import* 

[0040] from abaqusConstants import* 

[0041] from caeModules import* 

[0042] | executeOnCaeStartup () 

[0043] from abaqus import 为 定义 软件 为 ABAQUS 的 语句 ,from abaqusConstants 


import 语 句 含 义 为 导入 ABAQUS 中 的 
的 建 模 模块 ,execute0nCaeStartup 语 句 含义 为 对 建 模 模块 进行 初始 化 


[0044] 


身 管内 膀 结 构 如 图 


2 所 示 , 其 中 ,d 为 号 
坡 膛 长 度 ;d, 为 药 室 圆柱 部 直径 ,1, 为 药 室 圆 柱 部 长 度 ;ds 为 药 室 连接 锥 大 
连接 锥 长 度 ;d, 为 药 室 底部 直径 ,1 为 药 室 本 体 长 度 ;h 为 阴线 起 点 与 药 室 底 端 面 之 间 的 距 


常量 ,from caeModules import 语 句 含 义 为 导入 ABAQUS 


\ 一 /一 


运 们 。 


端 直径 ,1, 为 药 室 


fx 
= 
I1 


口径 ,L 8 


A ;kK 为 坡 膛 锥 度 。 此 外 将 膛线 条 数 记 为 n。 


[0045] | Egg 


HE ZEE AR 


3 所 示 , 其 中 ,d; 为 身 管 炮 口外 圆 直径 ,de 炮 口 制 退 器 安装 台阶 外 
管 尾 部 前 端 直径 ,d 为 身 管 尾部 外 圆 直径 ,Li 为 药 室 


圆 直径 ,d 为 身 管 中 部 外 圆 直径 ,d ,为 身 


底部 端面 与 炮 口外 


之 间 的 距离 ,L 为 身 管 尾部 长 度 ， 


出 退 器 安 厂 台阶 端面 之 间 的 距离 ,Le 为 药 室 底部 端面 与 吴 管 中 部 外 圆 端 面 


MN A L^ 
L 为 身 管 总 长 。 


[0046] ” 身 管 膛线 横 剖 面 图 如 图 4 所 示 , 其 中 ,a 为 阳线 宽 ,b 为 阴线 宽 ,t 为 膛线 深 。 
[0047] K= (d,-d) /1 
i 

[0048] hl *-L-lh-4l, 8 7 
[0049] F i Ze ih £ ET h £n PE 5 Pr s 185 Ht k RE 2k EH T 2 E ly TIE JT Bt ur, 组 
成 , 渐 速 段 曲线 方程 为 二 次 抛物 线 方 程 。 

MOM tan a, — tan ar, (= [tan a, | J [tan a (0s xx) 
[0050] 21 tanag,—tana,)  2(tano, -tanao) 

y,(x)- tana, :x ume me) LaL) 
[0051] ”其 中 ,y, G0 为 渐 速 段 曲线 方程 ,y, (x) 为 等 齐 段 曲线 方程 ,x 为 膛线 上 某 一 点 与 阴线 
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起 始点 之 间 的 距离 ,膛线 初始 缠 角 为 au, 炮 口 缠 角 为 ai ,膛线 总 长 为 L, 渐 速 段 长 度 为 1。 
[0052] 《将 au,ai,1,t,a 定 义 为 结构 参数 化 变量 ,并 赋予 初 值 。(ao\di 单 位 为 : 度 。1、t\a 单 位 
为 :mm) 

[0053] | a,—2.901,2,—7.778,1—4518,t—2.3,a—3.8 

[0054] ”将 身 管 内 膛 尺 寸 参数 中 的 d、d,、ds、dy、11、1,、1,、1l4sn 和 身 管 外 形 尺寸 参数 d,、d。、 
di dd、L LIL\L 均 设置 为 结构 参数 常量 。 
[0055] ”阴线 宽度 pbp 和 阴线 直径 di 与 结构 参数 化 变量 及 结构 参数 常量 之 间 存 在 以 下 关系 : 
[0056] ”阴线 宽度 :b= E -4，b/2 初 值 为 3.17, 膛 线条 数 n 王 48 ,阳线 半径 为 d/2 王 77.47; 
[0057] ”阴线 直径 :dj 二 2t+d,di/2 初 值 为 79.77。 

[0058] ”步骤 2 ,建立 映 管 及 弹丸 三 维 实体 模型 。 


[0059] ”根据 身 管 实际 结构 在 ABAQUS 软 件 中 按照 图 2 和 图 3 建立 直 简 身 管 三 维 实体 模型 。 如 
图 6 所 示 ,按照 图 4 在 建立 的 身 管 直 简 模型 炮 口 端面 绘制 膛线 截面 草图 ,膛线 条 数 为 48 条 , 膛 


线 分 界线 之 间 通 过 角度 进行 约束 ,角度 值 为 7.5 度 ;对 阳线 和 阴线 通过 半径 尺寸 进行 约束 ， 
明 线 半径 初始 值 为 79.77 ,阳线 半径 为 77.47; 炮 口 端面 草图 绘制 完成 后 ,采用 拉 伸 切除 命令 
得 到 身 管 直 膛线 轮廓 将 身 管 建 模 过 程 中 生成 的 fpy 中 相应 的 程序 段 提取 出 ,并 将 程序 段 中 
的 数值 采用 上 述 定 义 的 五 个 参数 化 变量 表示 ,得 到 直 膛 线 吴 管 实体 建 模 脚 本 程序 。 

[0060] ”弹丸 结构 示意 图 如 图 8 所 示 ,弹丸 为 底 四 弹 ,弹丸 由 弹 体 、 炸药、 引信 组 成 。 根 据 弹 
丸 的 实际 尺寸 建立 实体 模型 ,提取 出 建 模 过 程 中 rpy 文 件 中 相应 的 程序 段 ,将 程序 段 中 的 相 
关 数 值 采 用 上 述 五 个 参数 化 变量 表示 ,得 到 弹丸 实体 建 模 脚本 程序 。 

[0061] ”步骤 3 ,自动 划分 实体 模型 网 格 

[0062] ”采用 结构 化 网 格 划分 方法 对 直 膛 线 喘 管 实体 模型 进行 网 格 离散 ,建立 的 直 膛 线 母 
管 网 格 模型 如 图 7 所 示 。 弹 丸 有 限 元 网 格 模型 如 图 9 所 示 。 其 中 身 管 网 格 自动 划分 程序 如 下 
Bis 

[0063] ”设置 喘 管 网 格 种子 大 小 

[0064] | p—mdb.models['Model-1'].parts['SG'] 

[0065] | p.seedPart (size—5,deviationFactor-—0.1,minSizeFactor-—0.1) 

[0066] ”设置 网 格 单元 类 型 为 C3D8R 

[0067] elemType=mesh.ElemType (elemCode=C3D8R,elemLibrary=STANDARD, 
kinematicSplit-AVERAGE STRAIN,secondOrderAccuracy=OFF,hourglassControl= 
DEFAULT,distortionControl- DEFAULT) 

[0068] | c—p.cells 

[0069]  pickedRegions- (cells,) 

[0070]  p.setElementType (regions-—pickedRegions,elemTypes- (elemType,)) 

[071] ”步骤 4, 导 出身 管 有 限 元 模型 inp 文 件 ,Python 语言 读 取 单元 节点 并 按照 膛线 缠 角 
变化 规律 进行 偏 移 。 

[0072] “规定 炮 尾 指向 炮 口 方向 为 x 正 方向 , 竖 直 和 同上 为 y 正 方向 ,z 轴 正 向 由 右手 法 则 确 


c 


AC o 


[0073] 4.1) 定义 身 管 装配 体 , 在 装配 体 模块 定义 身 管 装配 体 模型 。 
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[0074] ^ a-mdb.models["Model-1'].rootAssemblytt&óE X 2x f I ER 

[0075] | a.DatumCsysByDefault (CARTESIAN) 楷 定义 身 管 坐标 系 为 全 局 坐标 系 

[0076] | p-—mdb.models['Model-1'].parts['SG' ] 检 定义 身 管 装 配 体 模型 为 Model-1 
[0077] a.Instance (name— ' SG-1' , part — p, dependent =0N) 

[0078] 4.2) 导出 直 膛 线 身 管 inp 文 件 , 将 身 管 装配 体 模型 Mode1-1 的 Job-SG.inp 文 件 导 
出 。 

[0079] mmdb.Job(name-— 'Job-SG' ,model— 'Model-1',description-—^",type— 
ANALYSIS,) 

[0080]  mdb. jobs[ Job-SG' ] .writeInput (consistencyChecking- OFF) 

[0081] 4.3) 计算 身 管 膛 线 缠 角 沿 炮 膛 轴线 的 变化 规律 , 渐 速 段 缠 角 正切 值 变 化 规律 为 
y O ;等 齐 段 缠 角 正 切 值 变 化 规律 为 y 和 , (x) : 


tan a, — tan a [tan a 
y (x)= o| x+ £ (0<x<1) 
I tan & — tan ar, 


[0082] 
y; (x) 2 tana (1 « x € L) 


[0083] ” 身 管 网 格 单元 某 一 节点 绕 炮 膛 转 过 的 角度 9 (x) 与 节点 距 阴 线 起 始点 之 间 的 距离 x 
存在 如 下 关系 : 


[0084] 0(x)= 


2x 2tanc 

— -y= eX 

nd d 

[0085] ”其 中 ,a 为 身 管 轴线 方向 某 一 单元 处 的 膛线 缠 角 ,n 为 膛线 的 缠 度 。 


a (a) 250). ,(o« y «j 


[0086] ”单元 节点 绕 炮 膛 轴线 偏 移 的 角度 : 


6, (x)= ae) x(Le x € L) 


[0087] ”其 中 ,X 为 吴 管 单元 节点 与 阴线 起 始点 之 间 的 距离 ,9 (x) ZEE Et Er T RE 
规律 ,9, (x) 为 等 齐 段 身 管 节点 偏 移 规律 。 

[0088] — 原 节 点 坐标 偏 移 一 定 角 度 后 的 新 节点 坐标 : 

[0089] 渐 速 段 身 管 节点 坐标 (y,z) 偏 移 9, 角度 后 的 节点 坐标 (y oz?) : 
jn = ycosÓ, — zsin, 


= ysin,  z cos0, 


[0090] ”其 中 ,8 为 吴 管 某 一 节点 在 渐 速 段 偏 移 的 角度 。 
[0091] 等 齐 段 身 管 节点 坐标 (y ,z) 偏 移 0 角度 后 的 节点 坐标 (y” ,2”) ， 
pe y cosÓ, —2zsing, 


ysinO, 4 z cos, 


[0092] ”其 中 ,9, 为 身 管 某 一 节点 在 等 齐 段 偏 移 的 角度 。 

[0093] 4.4) 按照 Haz 7l fg AS VER ETT Er ETT A 

[0094] ”如 图 10 所 示 , 线 膛 身 管 有 限 元 单元 节点 偏 移 过 程 为 : Python 程序 首先 从 炮 尾 中 心 
点 开始 读 取 身 管 单元 层 坐 标 和 节点 坐标 值 , 当 坐标 值 x, 小 于 阴线 起 点 与 药 室 底 端 面 之 间 的 
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距离 h 时 ,单元 节点 坐标 不 变 , 当 单元 层 x, 坐 标 值 大 于 h 且 小 于 h+1l 时 ,根据 膛线 缠 角 变化 规 
律 可 得 到 喘 管 渐 速 段 单元 节点 偏 移 后 的 坐标 (x0,y ^ z^) , 同 理 , 当 单元 层 x, 坐 标 值 大 于 h+1 
时 ,进行 等 齐 段 身 管 节点 的 偏 移 , 得 到 等 齐 段 节点 偏 移 后 的 坐标 (x0,y”,z”) , 随 着 单元 识别 
层 编号 i 的 增加 ,单元 节点 按照 缠 角 变化 规律 沿 着 身 管 轴线 方 品 依次 进行 偏 移 ,直到 单元 层 
x0 坐 标 值 大 于 喘 管 长 度 L. 时 ,号 管 单 元 节点 偏 移 结束 ,得 到 节点 侦 移 后 的 混合 膛线 叶 管 有 限 
元 模型 Job-SG_NEW.inp 文 件 , 其 中 Job-SG_NEW.inp 文 件 中 的 新 模型 为 Job-SG_NEW。 
[0095] 《进一步 地 ,作为 一 种 实施 方式 , 当 1=L 时 ,可 以 进行 渐 速 膛线 身 管 的 有 限 元 建 模 。 
[0096] ”进一步 地 ,作为 另外 一 种 实施 方式 , 当 1=0 时 ,可 以 进行 等 齐 膛 线 身 管 的 有 限 元 建 
模 。 

[0097] 又 5 ,将 单元 节点 偏 移 后 的 身 管 Job-SG_NEW. inp 文 件 通 过 Python 语言 读 取 并 妓 
i , 即 得 到 螺旋 膛线 映 管 有 限 元 网 格 模型 。 

[0098] 5.1) 导入 节点 偏 移 后 的 映 管 网 格 模型 ,相应 的 导入 程序 语句 如 下 所 示 : 

[0099] mdb.ModelFromInputFile (name= 'Job-SG NEW',inputFileName- 'Job-SG 
NEW. inp) 

[0100] ”mdb.ModelFromInputFile 目 的 在 于 读 取 模型 Job-SG_NEW 

[0101] 5.2) 导入 部 件 及 装配 体 属性 
[0102] import part## 导 入 部 件 属性 
[0103] import assembly# 导 入 装配 体 属性 

[0104] — 捧 将 模型 Job-SG_NEW 中 的 螺旋 身 管 部 件 PART-SG 导 入 到 模型 Model-1 中 

[0105] | mdb.models['Model-1'].Part ( PART-SG' ,mdb.models[ Job-SG NEW'].parts(' 
PART-SG']) 

[0106] — su ;E XL EFECREBOVEBETEPART-SG- 1 

[0107] | a-mdb.models[ Model-1'].rootAssembly 

[0108] | mdb.models['Model-1'].Instance ( PART-SG-1' ,mdb.models['Job-SG NEW'] 
.rootAssembly.instances['PART-SG-1']) 

[0109] ”节点 偏 移 后 的 身 管 有 限 元 模型 如 图 11 所 示 , 脱 线 结构 如 图 12 所 示 。 

[0110] ”步骤 6 ,定义 弹丸 及 身 管 的 耦合 关系 。 通 过 在 有 限 元 软件 ABAQUS 中 对 身 管 及 弹丸 模 
型 进行 材料 属性 定义 .装配 定义 、 分析 步 定 义 、 接 触 约束 定 义 、 载 荷 及 边界 条 件 加 载 操作 。 引 
信 、 炸 药 、 弹 带 分 别 与 弹丸 本 体 绑 定 约束 。 定 义 接触 属性 , 弹 带 、 前 定 心 部 与 映 管内 膀 之 间 设 
置 接 触 属性 , 切 回 摩擦 系数 取 0.1 ,接触 属性 采用 罚 函 数 法 描述 。 模 型 采用 0 号 装 药 , 弹 底 施 
加 由 内 弹道 方程 组 解 算出 的 压力 曲线 .边界 条 件 为 吴 管 尾部 全 约束 。 将 操作 过 程 中 软件 后 
台 上 自动 生 成 的 rpy 文 件 中 程序 语句 提取 出 来 ,得 到 喘 管 .弹丸 耦合 关系 施加 的 Python 程序 语 
句 。 

[0111] 其 中 创建 有 限 元 显 式 动态 分 析 步 的 相应 程序 段 如 下 所 示 : 

[0112] mdb.models[’Model-1’] .ExplicitDynamicsStep (name— 'Step-1' ,previous 一 ” 
Initial ,timePeriod=0.016, improvedDtMethod=ON) 

[0113] ”步骤 7 ,创建 静态 、 动 态 分 析 作 业 。 

[0114] ”复制 模型 并 创建 静态 分 析 步 

[0115] | mdb.Model (name =’ Model-2' ,objectToCopy=mdb.models[’Model-1°’]) 


© 


—— 


CN 112784458 A 说 HJ B 7/8 页 


[0116] — mdb.models['Model-2'].StaticStep (name— 'Step-1',previous-— 'Initial', 
maintainAttributes-True,nlgeom- ON) 

[0117] ”其 中 ,Model-1 为 动态 分 析 模 型 ,Model-2 为 静态 分 析 模 型 。 

[0118] — 7.1) 创建 静态 分 析 作 业 

[0119] — mdb.Job (name =’ Job-2' ,model — 'Mode1-2' ,description-—" , type — ANALYSIS, 
atTime-None,waitMinutes-0,waitHours-0,queue-- None memory — 90,memoryUnits-— 
PERCENTAGE,getMemoryFromAnalysis-True,explicitPrecision-SINGLE, 
nodalOutputPrecision-SINGLE,echoPrint-OFF,modelPrint-O0FF,contactPrint-— 
OFF,historyPrint-OFF,userSubroutine-—^,scratch-—" ,resultsFormat- -ODB, 


multiprocessingMode- DEFAULT , numCpus —4 , numDomains —4 , numGPUs — 0) 

[0120] 7.2) 创建 动态 分 析 作 业 

[0121] — mdb. jobs[ Job-2'] .waitForCompletion() 和 暂停 当前 脚本 执行 ,直到 当前 计算 任 
务 终止 

[0122] mdb.Job (name =’ Job-1' ,model= 'Model-1',description-" , type — ANALYSIS, 
atTime-None,waitMinutes-0,waitHours-0,queue--None memory — 90,memoryUnits-— 
PERCENTAGE ,explicitPrecision- DOUBLE PLUS PACK,nodalOutputPrecision-FULL, 
echoPrint-—O0FF,modelPrint-O0FF,contactPrint-—O0FF,historyPrint-OFF, 


userSubroutine-—" scratch" ,resultsFormat —0DB,parallelizationMethodExplicit 


—pOMAIN,numDomains-4,activateLoadBalancing-False,multiprocessingMode — 
DEFAULT , numCpus — 4) 

[0123] 7.3) 进行 重启 动 设置 ,设置 预定 义 场 ,实现 静态 分 析 作 业 与 动态 分 析 作 业 信 息 的 
交互 。 利 用 隐 式 求解 器 计算 重力 场 作用 下 的 初始 发 射 状态 , 设置 重启 动 分 析 ,将 重力 场 计算 
结果 作为 初始 条 件 施 加 在 动态 显 式 分 析 步 中 。 

[0124] ”静态 和 动态 分 析 步 之 间 设 置 预定 义 场 ,相关 程序 段 如 下 所 示 : 

[0125] | mdb.models['Model-1'].InitialState (updateReferenceConfiguration 一 ON， 
fileName =’ Model-1°,endStep=LAST_STEP,endIncrement =STEP_END,name =” 
Predefined Field-1',createStepName- 'Initial',instances-— instances) 

[0126] ”其 中 ,endStep 二 LAST_STEP 含 义 为 将 静态 分 析 步 的 最 后 一 个 分 析 步 计算 出 的 信息 
导入 到 动态 分 析 步 中 ,预定 义 场 被 命名 为 Predefined Field-1, 分 析 步 命名 为 Initial。 
[0127] 7.4) 保存 模型 

[0128] mdb.saveAs (pathName— ' 1") 

[0129] ”步骤 8 ,将 步骤 1- 步 又 7 中 的 Python 脚本 程序 进行 整合 后 ,得 到 弹 炮 耦合 参数 化 脚 
本 程序 。 采 用 ABAQUS 二 次 开发 技术 ,将 弹 炮 耦 合 有 限 元 建 模 过 程 通过 Python 脚本 程序 编写 
出 来 ,实现 弹 炮 耦合 有 限 元 模型 的 快速 建立 .建立 以 初始 缠 角 a 小 炮 口 缠 角 ai DR BEI ET. 
膛线 深 t 和 阳线 宽 a 为 变量 参数 的 弹 炮 耦 合 有 限 元 模型 时 ,如 图 13 所 示 , 首先 设 置 好 变量 的 
初始 值 ,然后 打开 ABAQUS 初 始 化 界面 ,点 开 RunScript 脚 本 运行 指令 ,如 图 14 所 示 ,选中 弹 炮 
耦合 参数 化 脚本 程序 py 并 点 击 运行 ,ABAQUS 软 件 将 按照 程序 代码 ,依次 完成 建立 实体 模型 、 
划分 网 格 、 定 义 并 赋予 材料 属性 、 设 置 接 触 与 约束 关系 、 施 加 载荷 与 边界 条 件 。 整 个 建 模 过 
程 全 部 交 由 脚本 程序 一 键 完 成 。 建 立 的 弹 炮 耦合 有 限 元 模型 如 图 15 所 示 。 
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[0130] 


以 某 大 口径 火炮 吴 管 内 膛 实 际 结构 参数 建立 的 弹 炮 耦合 有 限 元 模型 ,本 次 建 模 过 
程 所 花费 时 间 不 到 一 分 钟 , 与 传统 建 模 方法 相 比 ,运用 该 参数 化 脚本 程序 建立 弹 炮 耦合 模 
型 的 效率 会 极 大 提高 。 
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